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. Abstract 

' We obtain bi-Hamiltonian structure for a family of integrable systems on the sphere 

^ ^ , with an additional integral of third order in momenta. These results are applied to the 

O ' Goryachev system and Goryachev-Chaplygin top for which we give an explicit procedure 

^ ' to find the separated coordinates and the separated relations. 
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1 Introduction 



We address the problem of the separation of variables for the Hamilton-Jacobi equation within 
the theoretical scheme of bi-hamiltonian geometry. We want to learn to calculate bi-hamiltonian 
structure for a given integrable system on the Poisson manifold M with the Poisson bivector 
P and the Casimir functions Ca. The separation variables are naturally associated with this 
bi-hamiltonian structure of M itself [6] . 
^ ■ According to [HI [121 US] we will suppose that desired second Poisson bivector is the Lie 

f — I derivative of P along some unknown Liouville vector field X on M 

p' = cx{p) (1.1) 

o 



which has to satisfy to the equation 

[P',P'] = [Cx{P),Cx{P)]^0, ^ [Cj,{P),P]=0, (1.2) 

with respect to the Schouten bracket [.,.]. By definition (jl.ip bivector P' is compatible with a 
given bivector P, i.e. [P, P'] = 0. In geometry such bivector P' is said to be the 2-coboundary 
associated with the Liouville vector field X in the Poisson-Lichncrowicz cohomology defined by 
; P on M. 

^ ' From all the solutions X of the equation (|1.2p we have to choose partial solutions X such 

5^ ; that 

{H,,Hjy^O, z,j-l,...,n, (1.3) 

where Hj are the given integrals of motion and {., .}' is the Poisson bracket associated with the 
Poisson bivector P' (ll.ip . 

Obviously enough, in their full generality the system of equations (|1.2j|1.3p is too difhcult 
to be solved because it has infinitely many solutions labeled by different separated coordinates, 
see [13j . In order to get particular solutions we will use addition assumptions 

P'dCa = 0, a = l,...,k, (1.4) 

and some special ansatze for the Liouville vector field X. 

Equations (|1.4p say that we are looking for the Poisson bivector P', which has the same 
foliation by symplectic leaves as P. In general bivector P' could have some more Casimirs, so 
that their symplectic leaf could be smaller, but in the separation of variables method we have 
to consider equivalent foliations only [TTl [T^ [T3] . 
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2 A family of integrable systems on the sphere 

A Hamiltonian system is called natural if its Hamiltonian is the sum of a positive-definite kinetic 
energy and a potential. Natural Hamiltonian systems on cotangent bundles of closed surfaces 
admitting integrals polynomial in momenta are especially interesting 3 . In this Section we 
define some family of natural Hamiltonian systems on T*S^ with the quadratic and cubic in 
the momenta integrals of motion. 

Let two vectors J — {Ji, J2, J3) and x — {xi,X2,X3) are coordinates on the Euclidean 
algebra e*(3) with the Lie-Poisson bracket 

where Sijk is the totally skew-symmetric tensor. This bracket has two Casimir functions 

3 3 

Ci ^ \x\^ = Y,^l. C2 = (a:,J) = ^a;fcJfc. (2.2) 

k=l k=l 

Fixing their values one gets a generic symplectic leaf of e*(3) 

Oab- {x,J : Ci^ A, C2^B}, 

which is a four-dimensional symplectic manifold. 

At C2 = this symplectic manifold is equivalent to cotangent bundle T*S^ of the sphere 

5^ = {a;eR^ \x\^A}. 

For the Liouville integrability of the equations of motion it is enough to find one additional 
integral of motion iJ2, which is functionally independent of the Hamiltonian Hi and the Casimir 
functions. 

In this e*(3) coordinates the family of integrable systems on the sphere is defined by the 
following quadratic Hamiltonian 

Hi = J2 + Jl + {ia^ + /(X3)) Jl + m(x3)a;i -f 5(2:3) , (2.3) 

and cubic additional integral of motion 

H2 = -2aJ3(-aV| + Jl + Jl + /(xs) J3' + g{x3)) - n{x3)Ji - i{x3)xiJ3 . (2.4) 

Here 



(2.5) 



, . a^b n'jxs) , , njx^yjxs) 
5(2^3 = -7 — r^, m{x3) = , i[X3) = — r , 

2 3x3m(x3) - 2{A^ - xl)m'{x3) x^ejxs) ~ {A^ ~ 4y'{x3) 

f{x3) = 1 - 3a -a — ^ H — ^ , 

n{x3) n{x3) 

and 

n{x3) ^ -a{x3 + , a,5,ceR. (2.6) 
At C2 = integrals of motion Hi (|2.3p and H2 (|2.4|) are in the involution with respect to the 
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brackets in the following five cases: 

1. a — f3 — I, c = 0, n{x3) — —ax^, 



1 



1/3 



2. a — P——, c — 0, n(x3) = —axj 

3. Q; = /3=i, c = A, n(x3) = -a(2;3 + ("27) 

4. a = /3=i, ceM, n(a;3) = -a(x3 + c)^/^ , 

5. a = (3^^, c^A, n{x3) ^ -a{x3+ Ay/'^,. 
The corresponding Hamiltonian p.3p depends on parameters a, 6 and c: 
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rr(2) _ .2 , .2 , 4 2 aa;i & , . 

^1 - Ji + J2 + 0-^3 + -y^ + -TJz (2.8) 

O nit I rvt I 

^ - + + ^2 + 2(X3 + A) ) + {x^ + + (:,3 + ^)i/3 ' 

^ " '^i+'^^ + ^v X3+C 4(x3 + c)2;"^3 + (^^^^)i/2+a;3 + c' 

o-(5) _ 72 , t2 , /^H + ^ /2 , , ^ 

^ " ''^+''^ + ^6 8(X3 + A) j''3 + (x3+ A)3/4 + (^3 + ^)1/2- 

The Hamilton function ijj^'' describes the well-studied Goryachev-Chaplygin top IJ. In con- 
trast with this case the other systems have not at all physical meaning, however they may be 
interesting as some mathematical toys. The second integrable system with the Hamiltonian 
h'^^ was found by Goryachev [2J. The Hamilton function H^^^ was studied by Dullin and 
Matveev [1]. The third and fifth integrable systems with Hamiltonians iJ^-* and ijj^-* have 
been found in [5]. 

3 Construction of the Liouville vector field 

In this Section we describe how to obtain desired Liouville vector field X and the corresponding 
second Poisson bivector P' for our family of integrable systems. Remind, that smooth manifold 
M endowed with a pair of compatible Poisson bivectors P and P' is said to be loN manifold if 
one of the Poisson brackets is non degenerate [S] . 

Without loss of generality we can always put \x\ ^ A = 1 by using scaling transformation 
Xj — > A~^Xj and consider the unit sphere only. On its cotangent bundle T*S^ we introduce 
the following coordinates 

J1X2 - X1J2 

= arctan I — I, u = X3, = -J3, Pu ^ 2,2 — ' 

\X2 I 3Ti ~r 
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so that 



Ji = sin(0)p0 + Vl - cos((?!))p„, xi = \/l - 1*2 sin((?;)), 

V 1 — 



J2 = cos((/))p0 - Vl - sin(0)p„, ^2 = \/l - cos(0), 

In these coordinates u,(p,pu,P(f, initial Poisson tensor P associated with the brackets (|2.ip 
becomes canonical tensor on 2-dimensional symplectic manifold 



P = 



/ 1 \ 

1 

-10 

\ -10 0/ 



(3.1) 



which is non degenerate. 

Calculations in these coordinates are faster than calculations in the standard spherical 
coordinates (l>,0 — arccosw. Of course, using these coordinates on symplectic leaf of e*(3) we 
implicitly accepted assumptions (|1.4p for the initial Poisson manifold e*(3) with the Casimir 
functions ((2?2)) . 

The Hamiltonian (|2.3|) in these coordinates looks like 

Hi = (1 - u')pI + {ia" + f{u) + + - sin(0) m{u) + .g(«) 

= guPl + g^pl + V{u,(j,), (3.2) 



where 



(l-u^), g^ = 3a" + /(u 



u 
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1 -u^ 

According to [S] in order to solve equations (|1.2m.3p we will use polynomial in momenta 
ansdtze for the components of the Liouville vector field X — X^ dj : 

N k 

^' = EEyL.K'/')Pn->0- (3.3) 

fc=0 m=0 

For all the systems we put iV = 2, it means that X'^ will be generic second order polynomials 
in momenta Pu,P(p with coefficients y^„(u, 4>) depending on variables u and (j). 

Substituting this ansatz p.3p into the equations (|1.2j|1.3p and demanding that all the 
coefficients at powers of Pu and p^ vanish one gets the over determined system of 60 algebro- 
differential equations on the 24 functions J/j;m(u, (/>) which can be easily solved in the modern 
computer algebra systems. 

Below we discuss all the obtained solutions. 

3.1 Particular solution 

For the Goryachev-Chaplygin top [1\ with integral of motion h[^1 one get the following 
Proposition 1 At the first case a — (3 — 1 and n^x^) — ax^ there is one particular solution 
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depending on coordinate u only: 

/ M 



P' = 







(1-2m2)p^ 



(1 - u'^)pu 



1 



(3.4) 



y * * * 

The entries of P' are real functions on initial variables. 

By construction this Poisson bivector compatible with P (j3.ip and functions H12 are in bi- 
involution with respect to the corresponding Poisson brackets. So that the phase space T*S^ 
becomes semisimple toN manifold and the foliation defined by Hi 2 is separable in the Darboux- 
Nijenhuis variables [5]. 

The variables of separation qi^2 (the Darboux-Nijenhuis variables) are the eigenvalues of the 
recursion operator iV = P' P~^. They are simply roots of the following minimal characteristic 
polynomial of N 



AiX) =(dct {N-Xl) 



1/2 



= X"" - 2p4,X - (l ~ u^)pI 
= + 2./3A - J'^ - Jl - 



9 2 

^ pi 



(3.5) 



At & = these variables of separation have been found by Chaplygin in [T]. In initial e*(3) 
variables this bivector 



P' 



(0 


^3 


2^32^2 


-X2J1 


-X2J2 


X3J2 - 2X2 JsX 









XiJi 


X1J2 


2x1 J3 - 2:3 Jl 















—X1J2 + X2J1 











_t2 _ t2 
Jl J2 


-J3J2 


* 











J1J3 


V 










J 



(3.6) 



has been found in [T^ . It is easy to prove that P' has the same foliations by symplectic leaves 
as P at C2 = 0. 



3.2 Generic solution 



Using quadratic in momenta ansatz (j3.3p one got the generic solution depending on a pair of 
parameters fi and ly for all five cases of integrable systems 



Proposition 2 For integrable system i2.S^2.4^ equations il.S^l.S]) have the following solution 

fox y + AianP4, \ 

Z 4 ia^ptj, 

im(u)e'"^ 



P' = 



X 



V 



2^/T 




(3.7) 
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whe 



X 

y 



-2i/i(u + d) + 
X 



n{u)m{u) ' 
3a — 1 u 



l-v? 



P4> 



Here is a component of metric i3.2\} and 
• d = at cases 1,2,4, 



(l-u2) 



(1 -U2) 



• d — A at cases 3,5. 

The complex conjugated bivector P' is another solution of the same equations hl.2M.3^) . 

In contrast with the particular solution p.4p entries of this Poisson bivector are complex 
functions on initial variables and P'^^, depends on parameter a instead of parameter h. 

As above this Poisson bivector is compatible with P (13. ip and functions H12 are in bi- 
involution with respect to the corresponding Poisson brackets. So, the phase space T*S'^ be- 
comes semisimple uN manifold and the foliation defined by i?i.2 is separable in the Darboux- 
Nijenhuis variables, which are the eigenvalues of the recursion operator N = P'P^^ 6\. 

Summing up, we have found particular and generic solutions of the equations (jl.2j|1.3p 
for five integrable systems (|2.8p on the sphere with cubic integrals of motion. An important 
application of this result is the separation of variables for these systems. 



4 Separation of variables 

In this section we consider new separated variables and separated relations for the Goryachev 
system and for the Goryachev-Chaplygin top [I] in details. 



4.1 The Goryachev system 

In this case 

1 



fl - 1/3 « 

a = fi = n = -ax^' , m = 



2a 



3x 



2/3 ■ 



2/3 ' 



/=1, 



so that 



and 



rr(2) _ .2 . .2 , 4 2 , b 

til - ^1 + '^2 + 3 J3 + — + — 



H. 



(2) 



3^3 J? + j|+gJ3^ + ;;57^ 



1/3 J 2axi J3 
axr( Ji 



3x 



2/3 



If we put /i = and ly = — in P' (|3.7p one gets the following second bivector on T*S'^ 



I 



P' = 



^1/3 



1 /3 

iu^/'^p„ - 

1 — 



* 



3u2/3(l - u2)P" + 3(1-^2)2 





■P<p 



2wi/3VT^r^ 
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The variables of separation gi_2 (the Darboux-Nijenhuis variables) are eigenvalues of the recur- 
sion operator N = P'P~^, which are roots of the following polynomial 



^(A) = (A - q,)(X q^) = A2 + ^1/3 (^_^ _ ip^^ j X _ 
In the initial e* (3) variables the second Poisson brackets look like 



lae 



(4.1) 



{Xj , Xj } 
{Xl, Jl}' 

{X2.J2Y 
{Xl, J2}' 

{X2, J2}' 
{X3, Jl}' 



1/3 

T 4/3 T < 1/3 T 

X2J1 J2 4X3' X2J3 



30:3/^ xi+ix2 3(2:1 +1x2) 

• 4/3 T T .. 1/3 

iXg Ji a;iJ2 41X3 X1J3 



+ 12^2 3a;^^^ 3(a::i+ix2) 



i{x2 — 1x1X2 — 3a;§)J2 ^ 4ia;3^'^X2J3 
3a;3/^(a;i +1x2) 3(xi+ia:2) 
{xf + 1x1X2 — 3x3) Ji 4x3^'^a;i J3 



3x3'^(a;i +ix2) 3(xi+ia;2) 

xl^^{JlX2 - X1J2) f ixl^^{JiX2 - X1J2) 

{X3,J2} = -- 



Xl +1x2 

O'^ijkXk 
2X3^^ (xi + 1X2) 



X1+1X2 
i {J1X2 - xiJ2){Ji + i J2) 
3x3^^(2:1 + 1x2) 



The corresponding Poisson bivector P' is compatible with canonical Poisson bivector on e*(3) 
and satisfies to equations (|1.4p . so that P' has the same foliations by symplectic leaves as P at 
C2 =0. 

According to [6J the bi-involutivity of integrals of motion 



is equivalent to the existence of non-degenerate control matrix F such that 

2 



j(2) „{2)y 



P'dHP = Pj^P^jdH} 



i = 1,2. 



(4.2) 



In our case the control matrix F reads as 



/ X^/'(Ji+iJ2) , 2J3 



F : 



1 \ 



Q 2/3 
3X3' 



-ax 



2/3 



Xl -I- 1x2 

2xJ/'(Jl-HiJ2)J3 4J32 



\ 2(xi-fix2) 3(xi+ix2) 



9x 



2/3 



2/3 

-2J3 

3x^/^ / 



The Darboux-Nijenhuis variables qi^2 are simultaneously eigenvalues of the recursion operator 
and eigenvalues the control matrix. In our case they are the roots of the following polynomial 



^(A)=det(^^-AI) = A2 + "^^'^^^+^^^) 



Xl -I- 1x2 



A + 



2(xi -I- 1x2) 



The left eigenvectors of F, if suitable normalized, form the Stackel matrix S, which enters into 
a pair of the separated relations 



J2 S,,{q,,p,)Hf ~ U,{q,,p,) = , z = 1, 2. 



(4.3) 
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Here Ui are the Stackel potentials and pi,2 are variables conjugated to qi^2' 

Ui,?]} = Sij, {qi,?]}' = {qi,q]} = {qi,qj}' = {pi,Pj} = {pi,P]}' = O- 

Unfortunately, construction of the variables pi^2 is non-algorithmic procedure, which depends 
on the fortune and skilfulness In our case we can observe that 



qi + q2 



lU^'^Pu 



1-W2' 



and 



91 92 



2VT^ 
so that 

W,9i+92} = 0, {^0,9192} -i(?ig2, + 92} = iw^/^, {u,(7ig2} = 0. 

Integrating these equations with respect to p^{q,p) and u{q,p) we can easily get the following 
expressions for these functions 

3/2 

P0 = , and u 



qi - 92 

which yield the necessary definitions of the momenta 



-2i (gipi - 92P2) 



Pi 



where 



B{X) = i 



3(91 - 92) 
2 A 



3x 



2/3 



2 A 



Using the following relations for the second Poisson bracket 

{P4„ 9i + 92}' = 9i92, and {p^, 9152}' 

it is easy to prove that {qt^pj}' — Sijqi. 

In the separated variables the Stackel matrix S is equal to 



5* = 



1 

3i 



1 

3i 



1 

-1 



1 

-1 



291PI 2(72P2 



2 i a qipi 2 i a 52P2 



whereas integrals of motion have the generalized Stackel form at the second case ([2 



(2) 



/ (4pj + 6ip, -9)(3i-2p,)g: 



3ia^ 



3 16 



18p, 



iP] "^Pj 



i = 1,2. 



The corresponding separated relations (|4.3[) define two copies of the following algebraic curve 

^2 



4A' 



H = 2iaqiPi, A = . 



(4.4) 



In contrast with the usual Stackel systems it is non-hyperelliptic curve, which is trigonal alge- 
braic curve. It is interesting that one formally gets the real algebraic curve. But we have to 
bear firmly in mind that our separated variables 91.2 and pi^2 are complex functions on initial 
physical variables. 

If we put 1/ = and /i = 1 in P^^ l\3.7\i one gets another Poisson bivector on T*S^: 



I 



P' 



u 



* 



1 - 4u2 



3(1-^2)^ 



iM(4 - 1*2) 
3(1-U2)^"+ 3(1-^2)2' 





p<k 



1 a u 




2 

l/3pi0 
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In this case the Darboux-Nijenhuis coordinates are roots of the following polynomial 



A{X) = - iupu 



3(1 -u2)^ 



//3gi0 2iu 



2(2 + 



2vr 



3 W ' 9(1 „y2)^ 



In term of the previous separated variables this polynomial looks like 



2i 



2i 



-4(A) = A + — qipi A + — q2P2 



Now it is easy to prove that in the separated variables {q, p) bivector P' (|3.7p looks like 



P' 



-2 i a/i 



/ gipi \ 

q2P2 

-qiPi 

V -q2P2 / 



7^ 



( gi 

92 

-qi 

V -92 / 



(4.5) 



For the remaining four systems bivector P^^ (13. 7p has the similar form in the Darboux-Nijenhuis 
variables {q,p). 

4.2 The Goryachev-Chaplygin top 

Now we briefly discuss new separation of variables for the Goryachev-Chaplygin top [T] . Remind, 
that in this case the Hamiltonian reads as 



H 



(1) 



Ji + J2 + 4 J3 - axi + 



The new separated coordinates are roots of the polynomial 
AW 



\2 , 2 MP" , 3q!P0 ^ p^{uc+l) 

A + ft 



mn 77171 (u + c) 77177(7i + c)(m2 — 1) 



A + 



217111^^/1 — V? 



u 7/^(71^ — 1) 

whereas the conjugated momenta are equal to 



A + 



1 a e 



277VI - u2 



p,^B{\ = q,), where S(A) = i ( y - ^ 



These separated variables lie on two copies of the following algebraic curve 

,2 



C 



fi = 2ia qiPi, X ^ qi . 



(4.6) 



As for Goryachev model these separated variables 51.2 and pi^2 are complex functions on initial 
physical variables, which lie on the non-hyperelliptic algebraic curve. 

Of course, we can repeat similar calculations for the remaining systems and prove that 
they are related with different trigonal curves as well. As an example, for the DuUin-Matveev 
system [4j equation of motion are linearized on the following algebraic curve 



C : 



fi{fi'-X'-H[^^]=bX-H^'''-^, fi^2iaq,p„ X^q, 



r(4) 



4A' 



(4.7) 



if c = in the Hamihonian h['^^ (12^8)) . 

It will be interesting to get solutions of the equations of motion in term of the abelian 
functions for trigonal curves, as example see [5]. The other open question is construction of the 
Lax matrices associated with separated variables on trigonal curves. 
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4.3 The Jacobi method 



The Jacobi method consists of construction of the integrable system starting with some known 
separated variables and arbitrary separated relations. The method was originally formulated 
by Jacobi when he invented elliptic coordinates and successfully applied them to solve several 
important mechanical problems, such as the problem of geodesic motion on an ellipsoid, and the 
problem of planar motion in a force field of two attracting centers. In [7| Jacobi himself wrote: 
"The main difficulty in integrating a given differential equation lies in introducing convenient 
variables, which there is no rule for finding. Therefore, we must travel the reverse path and after 
finding some notable substitution, look for problems to which it can be successfully applied." 
In our case this notable substitution is described by the following 

Proposition 3 At any a, c and a transformation 

- c, p0 = , (4.8) 



gi - 92 / qi - qi 

21n(M + c)(3a- 1) +ln ' '^^i'^2(^ 



/ 3iQ i(MC+l) \ ., ^7r,-^ 

\u + c (w + c)(m^ — 1)/ 
is canonical transformation. 

Now in order to get some integrable system on T*S^ we could take two copies of any algebraic 
curve C defined by equation $(/i, A, Hi, H2) = and solve the corresponding separated relations 

^{Hi,Xi, Hi,H2) ^0, ^i.i^2iaqiPi, Xt ^ qi 

with respect to integrals of motion Hi and H2, which will be in bi- involution |13j . As above 
the main problem is that change of variables (14. 8|) is the transformation over the complex field 
C and if we want to get real functions i?i.2 on the initial variables on T*S^ we have to start 
with the very special algebraic curves C, for instance with (|4.4l4.6p or (|4.7p . 

Example: Let us consider the following deformation of the algebraic curve (14. 6[) : 



C: ^i^^^i' + ^iX-Hl)+p^^^ + ^l)ix + ^i) = bX-H2- — . (4.9) 

Solving the corresponding separated relations with respect to Hi 2 one gets integrals of motion 
in the bi- involution. After canonical transformation (|4.8p at c = and a — 1 we obtain the 
following complex Hamiltonian 

Hi = {l-u )p„ + — - a sm(0) V 1 - + — + pp^ — \ , 

1 — u 

which after an obvious additional shift of momenta (canonical transformation) 

Pu ^ Pu Pu r\ 

2u 

becomes the real Hamiltonian for the Goryachev-Chaplygin gyrostat 

Hi = h[^^ + pp^ ^ Jf + Jl + 4J2 - pjg -axi + \. 
More complicated deformations of the algebraic curve (|4.6p . such as 

2 

pi^p? + pX- Hi^ + cipX + C2y? + C3/i^ = bX - H2 ~ 

and the corresponding additional shifts of the momenta Pu lead to another generalizations of 
the Goryachev-Chaplygin top, which was obtained in [3]. 

Of course, we can get similar deformations for the remaining four systems in (|2.8p as well. 
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5 Conclusion 



We found two-parametric Poisson bivector l|3.7p . which is compatible with canonical Poisson 
bivector on cotangent bundle T*S^ of two-dimensional sphere. The quadratic and cubic inte- 
grals of motion (|2.3I2.4[) for the five integrable systems on the sphere are in bi-involution with 
respect to the corresponding Poisson brackets. 

The eigenvalues of the corresponding recursion operator are the separated coordinates. 
For the Goryachev system and Goryachev-Chaplygin top we give an explicit formulae for these 
separated variables and the corresponding separated relations. 

The research was partially supported by the RFBR grant. 



References 

[1] D.N. Goryachev, On a Motion of a Heavy Rigid Body About a Fixed Point in the Case of 
A = B ^4C, Mat. sbonik kruzhka lyub. mat. nauk, vol. 21, no. 3, pp. 431-438, 1900. 

S.A. Chaplygin, A New Partial Solution of the Problem of Rotation of a Heavy Rigid Body 
About a Fixed Point, Trudy otd. fiz. nauk Mosk. obshch. lyub. estest., v. 12, n.l, p. 1-4, 
1904. 

[2] D.N. Goryachev, New cases of integrability of Euler's dynamical equations, Warshav. Univ. 
Izv., V.3, p.1-15, 1916. 

[3] G. Darboux, Legons sur la Theorie Generale des Surfaces, 1898. 

[4] H.R. DuUin, V.S. Matveev, A new integrable system on the sphere. Mathematical Research 
Letters, v.ll, p. 715-722, 2004. 

[5] J.C. Eilbeck, V.Z. Enolski, S. Matsutani, Y. Onishi, E. Previato, Abelian functions for 
trigonal curves of genus three. International Mathematics Research Notices, v. 2007, Art. ID 
140 (38 pages), 2007. 

[6] G. Falqui, M. Pedroni, Separation of variables for bi-Hamiltonian systems. Math. Phys. 
Anal. Geom., v.6, p.139-179, 2003. 

[7] C.G.J. Jacobi, Vorlesungen iiber Dynamik, Georg Reimer, Berlin, 1866. Jacobi's lectures 
on dynamics given in Konigsberg 1842-1843 published by A. Clebsch. 

[8] A.V. Tsiganov, On a family of integrable systems on with a cubic integral of motion, 
J. Phys. A, Math. Gen. v.38, p.921-927, 2005. 

[9] A.V. Tsiganov, On natural bi-hamiltonian systems, Teor. Math. Phys., v. 149(2), p. 1437- 
1456, 2006. 

[10] A.V. Tsiganov, A family of the Poisson brackets compatible with the Sklyanin bracket, J. 
Phys. A, V.40, pp.4803-4816, 2007. 

[11] A.V. Tsiganov, Separation of variables for a pair of integrable systems on so* (4), Doklady 
Math., v.76(3), p.839-842, 2007. 

[12] A.V. Tsiganov, On bi-hamiltonian structure of some integrable systems on so* (4), J. 
Nonlinear Math. Phys., v.l5(2), p.171-185, 2008. 

[13] A.V. Tsiganov, On bi-hamiltonian geometry of the Lagrange top, J. Phys. A: Math. Theor., 
v.41, 315212 (12pp), 2008. 

[14] H.M. Yehia, On certain two-dimensional conservative mechanical systems with a cubic 
second integral, J. Phys. A: Math. Gen., v.35, pp. 9469-9487, 2002. 



11 



